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Abstract 

We perform the stochastic quantization of scalar QED based on a generalization of 
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1 Introduction 

The stochastic quantization scheme of Parisi and Wu [1] has been applied to QED since 
many years. Nice agreement with conventional calculations was found in several explicit 
examples (for reviews see e.g. [2, 3]), a general equivalence proof so far was lacking. 

The main idea of "stochastic quantization" is to view Euclidean quantum field theory 
as the equilibrium limit of a statistical system coupled to a thermal reservoir. This system 
evolves in a new additional time direction which is called stochastic time until it reaches 
the equilibrium limit for infinite stochastic time. In the equilibrium limit the stochastic 
averages become identical to ordinary Euclidean vacuum expectation values. 

There are two equivalent formulations of stochastic quantization: In one formulation 
all fields have an additional dependence on stochastic time. Their stochastic time evolu- 
tion is determined by a Langevin equation which has a drift term constructed from the 
gradient of the classical action of the system. The expectation values of observables are 
obtained by ensemble averages over the Wiener measure. 

Equivalently one has a Fokker Planck equation for the probability distribution char- 
acterizing the stochastic evolution of the system. Now expectation values of observables 
are defined in terms of functional integrals with respect to the stochastic time dependent 
Fokker-Planck probability distribution. The equilibrium limit of the probability distribu- 
tion provides the Euclidean path integral density. 

One of the most interesting aspects of this new quantization scheme lies in its rather 
unconventional treatment of gauge field theories, in specific of Yang-Mills theories. We 
recall that originally it was formulated by Parisi and Wu [1] without the introduction of 
gauge fixing terms and without the usual Faddeev-Popov ghost fields; later on a modified 
approach named stochastic gauge fixing was given by Zwanziger [4]; further generaliza- 
tions and a globally valid path integral were advocated in [5, 6, 7, 8]. 

The main difficulty for providing an equivalence proof in the case of QED appears 
to be a rather nontrivial topological obstruction; all previous attempts failed in the past 
years to identify the standard - gauge fixed- QED action as a Fokker-Planck equilibrium 
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distribution. 

In this paper we introduce new modifications of the original Parisi-Wu stochastic 
process of QED, yet keeping expectation values of gauge invariant observables unchanged: 
The modified stochastic process not only has damped flows along the vertical direction 
[4] but also is modeled on a specific manifold of gauge and matter fields with associated 
flat connection [5, 6, 7, 8]. It is precisely in this case that the standard - gauge fixed - 
QED path integral density can indeed be identified with the weak equilibrium limit of 
the underlying Fokker-Planck probability distribution. 

In section 2 the geometrical setting for QED is introduced and the associated bundle 
structure of the space of gauge potentials and matter fields is summarized. We introduce 
adapted coordinates, the corresponding vielbeins and metrics. 

The generalized stochastic process for QED is presented in section 3, section 4 is 
devoted to the derivation of the conventional QED path integral density as the equilibrium 
solution. 

2 The geometrical setting of QED 

In this section we present the major geometrical structures of QED. We collect in a 
somewhat formal style all the necessary ingredients which are needed for a compact and 
transparent formulation of the stochastic quantization scheme of QED. 

2.1 Gauge Fields 

Let P — > M be a principal [/(l)-bundle over the n-dimensional boundaryless connected, 
simply connected and compact Euclidean manifold M. The photon fields are regarded as 
elements of the affme space A of all connections on P. The gauge group Q is given by 
g = C°°(M;C/(1)) with Lie algebra LieQ = C°°(M;u(l)); here Lie (17(1)) = u(l) = iR. 
The action of Q on A is defined by 

A 9 = A + g' l dg. (2.1) 
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Let us define the subgroup Go C G where Go = Q/U{1) denotes the group of all gauge 
transformations reduced by the constant ones. Since Go acts freely on A we consider the 
principal Go bundle n : A — > M. = A/Q. One can prove that A — > M. is trivializable, a 
global section a : M. ^ A being given by <r([A]) = 1 . Here cj(A) G £ is defined by 

uj{A) = exp[&- l d*{A - A )} (2.2) 

where A denotes the invertible Laplacian which is acting on the Lie algebra LieGo ; 
A G A is a chosen fixed background connection. Note that u(A) fulfills 

u(A°)=u(A)g (2.3) 

with g G Go- 

2.2 Matter Fields 

In order to discuss scalar matter fields we chose a representation p of g £ Go on 
the vector space V = C ; p(g) simply denotes multiplication with g. We consider the 
associated vector bundles E = P x p V on M. Scalar fields are described by appropriately 
chosen sections of E. In the following we denote by T the space of scalar fields. 
The action of Go on := (A, <fi) e A x J 7 is given by 

& = (A 0) _^ = (A 9 , = (A + fT 1 ^, <T V). (2.4) 

We remind that A x T — .4 Xg JF is trivializable iff ^4 — ^ .M is trivializable. Indeed, 
using the previous construction of u>(A) we obtain a global section a : Axg J 7 — > Ax F 
by defining 

<j{[A,<t>\) = (a([A]),r {Arl ) = (A^)- 1 ,^)- 1 ). (2.5) 

The tangent space of the configuration space .4. x T is given by T^^.A x T) = 
fl 1 (M; i R) x JF. Noting that G JF by construction transforms equivariantly we ob- 
tain a (/-invariant Riemannian metric h : T(A x JF) x T(A x JF) — > R by defining 

V,,) ((7>J), (t>J)) = (r\,rl) + ( V J, V J>. (2.6) 
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Here 

(a, (3) = - [ (a A *(3 + a A */3) (2.7) 
2 Jm 

and a,f3 G TVl, or TjF, respectively ; * is the hodge operator on M, ct denotes complex 
conjugation of a. 

2.3 Adapted Coordinates 

Let the globally defined gauge fixing surface £ in A x T be defined by 

Y 1 = ima = {(B,^) eAx F\(B,4>) = (A^' 1 ,^' 1 )} (2.8) 

where u is given by (2.2). Note that B and ip are invariant under the action of Q which 
trivially follows from (2.4) and from (2.3); B satisfies the "gauge fixing condition" 

d*(B - Aq) = 0. (2.9) 

We define the adapted coordinates ^ = {B, ip, g} via the bundle maps X '■ £ x Q — > 
Ax J 7 and x^ 1 '■ A x T — > £ x £? , where 

X {B^g) = {B\^) (2.10) 

and 

X- 1 ^^) = (a(L4,0]),c(A)) = (A^)-\r (A) ~\^))- (2.H) 
The differentials Tx and Tx" 1 are calculated straightforwardly (compare also with [7]) 

T X (( B ,v^Y g ) = (( B + d9 g (Y g ),g- 1 (v^-9 g (Y g )<S>)), (2.12) 

as well as 

T X -\T A ,Vt) = (Pr A ,u(A)(v (l> + (A- 1 d*r A ) ( f ) ),u(A)A- 1 d*r A ). (2.13) 

Here (Cs,ity) G TgE, F 9 G T 9 £ and (ta, ity) G T(a,</,)(^4 x JF). The Maurer-Cartan form 
on Q is denoted by 9, P is the transversal projector P = 1 — dA~ l d*. From the differ- 
entials Tx and Tx^ 1 we read off the vielbeins e\ = and their inverses -E 1 ^ = ||£ 
corresponding to the change of variables = {B,ip,g} <-> = {A, 0}. 
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Above we defined a Riemannian structure on Ax T in terms of the Go invariant metric 
h. In adapted coordinates this metric is given now as the pullback G = x*h] explicitly 
we obtain 

= (Cb^Q + (e g (Y}), (a + M 2 )e 9 (Y*)) + („},„}> - (e 9 (Y g ^,vl) - (v^e,^). 

(2.14) 

Here (£#, v^), v^) G T( Sj ^,)E and Y^ 1 , F fl 2 G T 9 £. In matrix notation we have G = e 
and G^ 1 = ££t. The determinant of G is given by det G = A . 

3 Parisi-Wu Stochastic Quantization 

For scalar QED we have 

S inv = (D A (p, D A ip) + m 2 ((/9, ip) + (3.1) 

where D A f — (d — A)ip and F = dA. The Parisi-Wu Langevin equations are given by 

dA= JJ^ ds + du (3>2) 
5.4. 

^ = JJ^L ds + rfy (3_ 3 ) 

0<£> 

where the Wiener increments fulfill 

dU dU = 2ds, dV dV = 2ds. (3.4) 

This can be summarized by 

d$ = -?*^ds + d£, d£dg = 2'lds. (3.5) 

Using Ito calculus [9, 10] we transform the Parisi-Wu Langevin equations into adapted 
coordinates 



d^ = 



_! 5Sj nv 5G 1 
5^ 5^ 



ds + d(, (3.6) 
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where 

d(d^ = 2 G~ l ds. (3.7) 

The use of adapted coordinates allows to disentangle the complicated dynamics of gauge 
independent and gauge dependent degrees of freedom; it will be of great value later on. 
For completeness we note the Fokker-Planck operator in the original variables 



8 Si, 



5$ 



_6_ 
+ 5$ 



as well as in the adapted coordinates 



6 Si, 



5^ 



+ 5$> 



(3.8) 



(3.9) 



We remark that in the case of the Parisi-Wu processes diffusion along the vertical di- 
rection takes place and no equilibrium distribution is approached. Thus a Fokker-Planck 
formulation of the Parisi-Wu stochastic quantization scheme is impossible: The gauge 
invariance of the action Si nv is leading to divergencies along the vertical directions when 
trying to normalize the Fokker Planck density. 



4 Generalized Stochastic Quantization 

4.1 Geometric Obstruction 

Our equivalence proof relies on specific allowed modifications of the metric on the field 
space, which governs the stochastic process. These modifications correspondingly are 
implying changes of the associated Fokker-Planck operator. We are going to show that this 
can be achieved in such a way that the resulting Fokker-Planck operator has a positive 
kernel and is annihilated on its right by the standard gauge fixed QED path integral 
density. In order for this to be the case, however, a certain integrability condition for the 
drift term of the considered stochastic process has to be fulfilled. Surprisingly similar as 
in the pure Yang-Mills theory also in the abelian QED case there appears a violation 
of this condition; it is only after a nontrivial modification of the underlying stochastic 
processes (see next subsection) that this obstruction can be overcome. 
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Proceeding step by step we first note (see Zwanziger [4]) that a damping force along 
the gauge orbit has to be introduced in order to maintain the probabilistic interpretation 
of the Fokker-Planck formulation. Although knowing that this additional force will not 
alter expectation values of gauge invariant quantities it is disappointing to observe that 
due to its presence the standard - gauge fixed - QED action will never annihilate the 
Fokker-Planck operator on its right side due to the following reason: 

We recall that the bundle metric h(A,<j>) on the associated fiber bundle AxT — > AxqT 
which is invariant under the corresponding group action gives rise to a natural connection 
7, whose horizontal subbundle Ti is orthogonal to the corresponding group. The horizontal 
bundle H[A x T\ 7] with respect to 7 is defined by H[A x T\ 7] _U V(A x JF), where the 
orthogonality is with respect to h(A,<t>)- Elements of the vertical bundle V(A, T) — > AxgT 
are given in the form 

Z^AA) = |u(A e ^0 e ^) = (df,-^), (4.1) 
where £ G C°°(M;iR). The orthogonal span with respect to the vertical bundle fulfills 

d*r A + 7^00 _ v 4> ( t ) ) = °> t a^ T a A, Vtf, G T^T (4.2) 

which follows from 

= h {A ,4>) ((^4,^),Z f (A,0)) = h( A ,4>) ((TA) v <t>) > (<%,-&)) = {TA,d£} - {v^^4>). (4.3) 

Explicitly we can prove that 7(^,0) (TAj v^) G C°°(P;iR) 

1(a,4>){ta, V+) = (A + r A + ^0 - v<t> 4>)] (4.4) 

defines a connection induced by h(A,<t>) in the principal bundle Ax T — > -4xgjF and is {/ (1) 
invariant. Calculating its curvature f2 ((r 1 , v 1 ), (r 2 , t> 2 )) we find that it is nonvanishing and 
given by 

A = (A + |0| 2 )- 1 (t; 2 + t; 2 0)(A + |0| 2 )- 1 [rfVA + i(t;V-^ 1 0)] 
-(A+|0| 2 )- 1 (t; 1 + t;V)(A + |0| 2 )- 1 

+ [d*r 2 A + - v 2 0)}(A + l^l 2 )' 1 ^ V - v\ 2 ) (4.5) 
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As a consequence [5, 6, 7, 8] there does not exist (even locally) a manifold whose tangent 
bundle is isomorphic to this horizontal subbundle. Specifically this implies that any vector 
field along the gauge group cannot be written as a gradient with respect to the metric 
h(A,<t>)- The total drift term - containing the extra vertical force term - thus can never 
arise as derivative of the standard gauge fixed QED action; the Fokker-Planck operator 
can never be annihilated on its right by the standard QED path integral density; an 
equivalence proof presently cannot be given. 

4.2 The Induced Field Metric with Flat Connection 

The crucial observation in [5, 6, 7, 8] is to consider a larger class of modified stochastic 
processes than considered so far, yet always keeping expectation values of gauge invariant 
observables unchanged: One introduces not only the extra vertical drift terms as discussed 
above but one also modifies the Wiener increments by specific extra terms and introduces 
extra so called Ito-terms, correspondingly. 

The idea is to view the new terms multiplying the Wiener increments as vielbeins 
giving rise to the inverse of a yet not specified metric on the space A x T . The appearance 
of this metric induces a specific connection with a potentially analogous obstruction as 
discussed above. A necessary requirement to overcome this obstruction is therefore that 
the corresponding curvature has to vanish. The question how to find such a metric is 
reduced to the question how to find a flat connection. 

Indeed, there exists a flat connection 7 in our bundle. This connection is the pull-back 
of the Maurer-Cartan form 6 on Q via the global trivialization and prg 

7 = (X" 1 >5*)(A0 = A- I cr, (4.6) 

where prg is the projector S x Q — > Q. The projector onto the horizontal subbundle 
H[A x JF; 7] with respect to 7 is given by 

P = 1 - D A % (4.7) 

We see that the horizontal subbundle H is orthogonal to the gauge orbits with respect 
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to the induced field metric; in particular the gauge fixing surface is then orthogonal to 
the gauge orbits. 

In the adapted coordinates the induced field metric is denoted by G = e 1 e . The just 
discussed orthogonality condition of the gauge fixing surface and the gauge orbit with 
respect to the induced field metric is transformed into simply 

(GT 1 ) 130 = (G _1 ) eE = 0, where G' X =E& with Ee = l. (4.8) 



This condition is fulfilled provided E is defined as 

( E*\ 



E = 



(4.9) 



V e b J 

To complete our discussion we also have to specify the vertical drift term; it is related to 
the gradient of Sg, where we chose 



S g [g] = ^(d*g*e u ^,d*g*e u ^), 



(4.10) 



where 9 U ^ is the Maurer Cartan form on U(l). Note that in the original variables we 
obtain the standard background-gauge gauge fixing term 



( X - l *pr* g Sg)(A,<j ) ) = S g (uj(A)) = -(d*(A - A ), d*(A - A )). 



Summarizing we have 



where 



d^ 



^ SS tot 5G_ 
5^ 5^ 



ds + ( 

Stot = S inv + S g , and d(dC} = 2G~ 1 ds. 



(4.11) 
(4.12) 
(4.13) 



4.3 The Equivalence Proof 

It is easy now to prove for QED the equivalence of the stochastic quantization scheme 
with the path integral quantization. For the formulation in terms of the adapted coordi- 
nates \& = {B,tp,g} the associated Fokker-Planck equation is derived in straightforward 
manner 



ds 



= L[*]p[%s], 



(4.14) 
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where the Fokker-Planck operator L[ty] is appearing in just factorized form 

(4.15) 



5s tot m + _s_ 



5^ 5^ 



Due to the positivity of G the fluctuation dissipation theorem applies and the equilibrium 
Fokker-Planck distribution p cq [^l/] obtains by direct inspection as 

e -S to t[B,i>,g\ 



p eq m = -7— 

p -Si nv [B,ip] P -Sg[g] 

(4.16) 



^VBV^e~ s ^ B ^ Jg Vge- s ei9}' 

This result is completely equivalent to the standard background-gauge fixed QED path 
integral prescription. The additional finite contributions of the gauge degrees of freedom 
always cancel out when evaluated on gauge invariant observables. 

Similarly, in terms of the original variables $ = {A, 0} the Fokker-Planck equilibrium 
distribution is given by the standard background-gauge fixed path integral density 

e -Stot[A,4>] e -S inv [A,^]-Sg(w(A)) 

P ^ = J A VAV(pe- St ^ A ^ = J A VAV(pe- s ^ A ^- s s(^ A ))' 
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